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We present classical and quantum mechanical multiscale descriptions of lattice dynamics, from 
the atomic to the continuum scale, using atomic scale symmetry modes and their constraint equa- 
tions. This approach is demonstrated for a one-dimensional chain and a two-dimensional square 
lattice with a monatomic basis. For the classical description, we find that rigid modes, in addition 
to the distortional modes found before, are necessary to describe the kinetic energy, and obtain con- 
straint equations among these modes. Lagrangian equations, modified with the Lagrange multiplier 
terms, are solved for phonon dispersion relations without using displacement variables explicitly. 
The long wavelength limit of the kinetic energy terms expressed in terms of atomic scale modes is 
shown to be consistent with the continuum theory, and the leading order corrections are obtained. 
We also analyze the phonon in terms of symmetry modes, and find how the contribution of differ- 
ent symmetry modes varies depending on the phonon branch and wavevector. For the quantum 
mechanical description, we find conjugate momenta for the atomic scale symmetry modes. In di- 
rect space, graphical rules for their commutation relations are obtained. Commutation relations in 
the reciprocal space are also calculated. We emphasize that the approach based on atomic scale 
symmetry modes could be useful for the description of multiscale lattice dynamics, materials with 
electron-phonon coupling, and the dynamics of structural phase transition, which can be probed by 
time-resolved x-ray diffraction. 

PACS numbers: 61.50.Ah, 63.20.D-, 63.22.-m, 62.20.-x 



I. INTRODUCTION 

The dynamics at nanometer length scale has been 
a focus of recent attention.— In particular, materials 
with competing ground states, such as high tempera- 
ture superconducting cuprates and colossal magnetore- 
sistive manganites,^'^''* often show dynamic nanometer 
scale features, for example, stripes in cuprates^i^ and 
anisotropic correlations in manganitesjii^ Furthermore, 
recent advances in time-resolved x-ray technique have al- 
lowed experimenters to directly probe lattice dynamics 
in atomic scale. ^ It is believed that understanding these 
nanoscale features and their dynamics is essential to ex- 
plain macroscopic properties of these materials. 

For the description of mesoscopic scale domain struc- 
tures and their dynamics, phenomenological Ginzburg- 
Landau formalism has been very successful i^^*^ One of 
the keys for such a success is the use of symmetry in 
the definition of variables, which makes the selection of 
free energy terms self-evident. Motivated by the suc- 
cess of the Ginzburg-Landau approach for the continuum, 
symmetry-based atomic scale description of lattice dis- 
tortions has been recently proposed, and demonstrated 
for a two-dimensional square latticcji^ In this approach, 
atomic scale symmetry modes are defined on a plaquette 
of atoms, and are used to express potential energy terms 
associated with lattice distortions. This method has been 
used to understand atomic scale structure of twin bound- 
arieSjiS, antiphase boundaries and their electronic tex- 
tures,— strain-induced metal-insulator phase coexistence 
in manganiteS)^ superconducting order parameter tex- 
tures around structural defects^ and the coupling be- 



tween electronic nematic order parameter and structural 
domains in metamagnets near a quantum critical point 
Thus far, this approach has been used for static lattices, 
or the relaxation of lattice distortions introduced through 
the Euler methodji^ which does not require kinetic en- 
ergy terms. In the current paper, we present our study on 
how the approach based on atomic scale symmetry modes 
can be extended to describe lattice dynamics, within the 
scope of both classical and quantum mechanical formal- 
ism In Section II, we discuss how to express kinetic en- 
ergy term in symmetry modes, present our study within 
the formalism of classical mechanics, compare our result 
with the continuum results, and analyze the phonon 
spectrum in terms of symmetry modes. We formulate 
quantum mechanics in terms of atomic scale symmetry 
modes in Section III, with a summary given in Section 
IV. 



II. CLASSICAL FORMALISM 

A. One-dimensional lattice with a monatomic basis 

Using a one-dimensional lattice with a monatomic ba- 
sis shown in Fig. [1] we demonstrate mode-based descrip- 
tion of lattice dynamics. The displacements of atoms are 
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FIG. 1; One-dimensional lattice with a monatomic basis. 
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represented by Ui, where i being an index for the site. To 
be specific, we assume that the interaction between the 
nearest neighbor atoms are described by a spring with a 
spring constant K while other potential energy terms are 
negligible, as represented by the following Lagrangian, 

lK{u,+i-u,f, (1) 



where M is the mass of the atom. We take a two-atom 
unit as a motif for this lattice, and define the symmetry 
modes, e(i) and t{i), 

--{Ui+l - Ui), 



e(^) 



V2' 



1 



(2) 
(3) 



where a normalization factor is chosen according to the 
number of displacement variables in the definition. These 
modes are also shown in Fig. [2l 



FIG. 2: Symmetry modes for the one-dimensional chain in 

Fig.m 



and a well-known dispersion relation for the one- 
dimensional chain, 1^ 



(1 — cos k). 



(10) 



This result shows that the lattice dynamics can be stud- 
ied within the framework of atomic scale symmetry 
modes and their constraint equations, without using the 
displacement variables explicitly. 

Anharmonicity of one-dimensional chains is impor- 
tant to understand non-linear excitations, such as soli- 
tons, kink-solitons, intrinsically localized modes, and 
breathers li^ii^ Atomic scale modes, e and t, found here 
can be used to incorporate such anharmonicity into the 
Hamiltonian, which, along with their constraint equa- 
tions, would provide a formalism to study the dynam- 
ics of non-linear excitations in one-dimensional chains. 
In the next subsection, we demonstrate how the mode- 
based approach is applied to lattice dynamics for a two- 
dimensional square lattice with a monatomic basis. 



The two variables, e and t, correspond to the distortion 
and rigid translation of the motif, respectively. Since the 
two modes are defined from one physically independent 
displacement variable at each site i, these modes are re- 
lated through one constraint equation shown below in the 
reciprocal space and direct space, respectively. 



f{k) = {e'^ + l)e(fc) - {e'" - l)i(fc) = 
e{i -f 1) + e{i) - t{i + 1) + t{i) = 0. 



(4) 
(5) 



In terms of these modes, the Lagrangian in Eq. IH) is 
expressed in the following way 

ichain - E ['^(*)' + - (6) 



Two-dimensional square lattice with a 
monatomic basis 



Symmetry-based atomic scale description of lattice 
distortions for a two-dimensional square lattice with a 
monatomic basis, shown in Fig. [31 has been studied in 
Ref. 13, where dilatational ei, shear 62, and deviatoric 
63 modes, and short wavelength modes, and Sy, are 
defined, as shown in Fig. [H In terms of displacement 
variables u5 and uH, and and iy representing site in- 
dices, these distortional symmetry modes are expressed 



The result shows that introduction of atomic scale rigid 
modes, such as t, which are not considered in Ref. Il2l . 
allows kinetic energy term being expressed in a quadratic 
form. To obtain equations of motion, we modify Lchain 
with a Lagrange multiplier A(fc), as shown below. 



chain 



[e(fc)e(-fc)+i(fc)i(-fc)] 



ur 



i 



FIG. 3: Two-dimensional square lattice with a monatomic 
basis. 



--(2X)e(fc)e(-fc) + A(fc)/(-fc). (7) 

The Lagrangian formalism of dynamics leads to the two 
equations of motion. 



M 



e{k) + 2Ke{k)~\{k){e 



M 



t{k) + \(k){e 



-ik 



-ik 



1) 
1) 



0, 
0, 



(8) 
(9) 
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FIG. 4: Distortional symmetry modes of the motif for the two- 
dimensional square lattice with a monatomic basis in Fig. [31 
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as follows. 

ei{i) 



2^/2 



4 i 4+10 i+10 



e2& = 



-"^oi+"loi+"^ii+"lii 
1 



(11) 
(12) 

(13) 
(14) 
(15) 



Instead of Sx and Sy modes, the following s+ and s_ 
modes can be also used. 



2\/2 ^ * * 

+u5 - uH + u? + 

i+Ol j+Ol i+ll j+ll 

e3(?) = -^(-uS + uy.+uS +ul 
2\/2 ^ * * * 

-u?, „, - + M?, , - 

J+Ol i+Ol i+ll i+ll 



S2;(?) = - fu5 — u? 
^ ^ 2 V « 



'i+10 


- "toi - 




-11 


T+io 


i+Ol 


24 


-11 



s+{i) = -^[s^(i) -f Sy(i)], 



(16) 
(17) 



These five modes have been used to express various har- 
monic and anharmonic potential energy terms but 
are not sufficient to represent kinetic energy terms in a 
simple quadratic form. For example, an expression has 
been converted back to the displacement variables to ob- 
tain the phonon spectrum shown in Fig. 2 in Ref. [T^ . 

In current work, we show that additional modes asso- 
ciated with the rigid motion of the motif, similar to the 
mode t in the one-dimensional chain, allow a formalism 
entirely based on symmetry modes without resorting to 
displacement variables. Three rigid modes for the two- 
dimensional square lattice are shown in Fig. [S] and are 
defined as follows. 

^^(^ - ^("?+-^io+"l,oi+-^ii)' (18) 

- |("f + "lio+-loi+-lii)' (19) 



^ ' 2^2 V * 



I i+lO j + 10 



i+Ol i+01 2+11 i+ll 



(20) 



The first two modes, tx and ty, correspond to rigid trans- 
lations of the motif along x and y direction, and the mode 
r represents a rigid rotation of the motif. Following i+ 
and t_ modes can be also used as alternatives to tx and 
til • 



1 



= -^[tx{i) + ty{i)i 

t-(?) = -^[tx{l)~ty{l)]. 



(21) 
(22) 



FIG. 5: Rigid symmetry modes of the motif for the two- 
dimensional square lattice with a monatomic basis in Fig. |31 



Straight-forward expansion shows that the kinetic en- 
ergy of the lattice is expressed in terms of the eight sym- 
metry modes in the following quadratic form, with M 
being the mass of the atom. 



(23) 



+ Sy(if +ix(if +iy{if +f(lf 



(24) 



As discussed in Ref. 12:, constraint equations are found 
from the relations between symmetry modes and dis- 
placement variables in the reciprocal space. We first rep- 
resent (uS, m|) in terms of (sa;(fc), Sy{k)) by inverting 
the linear relations between them, and replace in the ex- 
pressions with other modes, which lead to six constraint 
equations 
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sin ^ cos ^Sx{k) + cos ^ sin ^Sy{k) - 

-V2isin^sin^ei(fc) = 0, (25) 
cos ^ sin ^Sx{k) + sin ^ cos ^Sy{k) — 

-\/2isin^sin^e2(fc) = 0, (26) 
sin !f cos !fsx{k) - cos !f sin ^Sy{k) - 

-^/2^sin^sin^e3(fc) = 0, (27) 
cos ^ sin ^Sx{k) — sin ^ cos ^Sy{k) + 

+V2i sin ^ sin ^r{k) ^0, (28) 



2 2 

cos%^cos^S:r(fc) +sin^sin^t:r(fc) = 0, (29) 
cos ^ cos ^Sy{k) + sin !f sin ^ty{k) = 0. (30) 



Modified Lagrangian for the square lattice is now 

6 



(31) 



where V is the potential energy, A„ (k) are Lagrange mul- 
tipliers, and fn{k) = are the six constraint equations, 
Eqs. (^5)1 - ([50)1 . By solving the Lagrangian equations, 
we find dynamic properties of the lattice in terms of 
the symmetry modes. As with the Ginzburg-Landau 
approach being useful for the description of mesoscale 
dynamics, we expect that the approach based on atomic 
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scale symmetry modes would be useful for the description 
of atomic scale dynamics, particularly, when anharmonic- 
ity plays an essential role. In the next two subsections, 
we present two applications of the formalism developed 
in this subsection. 



C. Comparison with continuum description of 
lattice dynamics 

We compare the atomic scale theory developed in the 
previous subsection with an existing continuum theory 
of lattice dynamics. Either by using the definitions, 
Eqs. ifTTj) and p3|) . or by using the constraint equa- 
tions, we express the kinetic energy for the square lattice, 
Eqs. ([23]) and (f24|) . in terms of ei and 63, 

t s— 1,3 s'— 1,3 



to in the continuum theory. 



ei{i) 


1 


+ V,uH], 


(37) 






+ V^u5], 


(38) 


eaii) 






(39) 


r{i) 






(40) 






? 


(41) 


Sy{i) 
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(42) 


tx{i) 


= 2m5, 




(43) 


ty{i) 


= 2mH. 




(44) 



In fc ^ limit, the correspondence of these modes to the 
displacements u are 



where 



711 (fc) = 733 (fc) 
7l3(fc) = 731 (^) 



1 — cos kx COS ky 

sin^ kx sin^ ky 
COS kx — cos ky 
sin^ kx sin^ ky 



(33) 
(34) 



^7 

61,62, 63, r ku, 



(45) 



The comparison shows that our approach is a natural ex- 
tension of the continuum theory to the atomic scale, and 
is suitable for multiscale description of lattice dynamics. 



To compare with a continuum theory, we take the long 
wavelength limit, and obtain the following leading order 
term for 755/ , 



D. Phonon analysis using symmetry modes 



ky ky 



91,2 h.2 



■il^x'^y 



k^ k^ ~\- k,. 



\ 



91,2 h.2 



■^'^x'^y 



\ 



J 



(35) 



We analyze phonon modes in terms of atomic scale 
symmetry modes for the square lattice with a harmonic 
potential shown below. 

= E j^iei(^)' + \A2e2{lf + ^^363(1')' 



This term is identical to Eq. (3.12a) in Ref. MT* (63 here 
corresponds to 62 in Ref. ilL). which Lookman et al. have 
used as continuum kinetic energy to study underdamped 
dynamics of strains in proper ferroelastic materials. The 
next order term to the above continuum limit is as fol- 
lows. 



V 



1 

12 



k^ ■ 

Ky 



1,21,2 
^x'^y 



\ 



Sk'^ky 



1 

12 



1,2 1.2 

^x'^y 



(36) 



This term, or better Eqs. ((33)) and ((34|) . can be used to 
study the dynamics of proper ferroelastic materials on the 
atomic scale. The following long wavelength limit of the 
atomic scale modes shows directly what they correspond 



+ -B[Sx{lf + Sy{lY 



(46) 



By solving the Lagrangian equations Eq. (|3ip . we find 
dispersion relations without using displacement variables 
explicitly, 

Mw^ = B(l - cosfc3;)(l - COSfcy) 

+ — {Ai + A2 + — cos kx COS ky) 



± 



-(^1 - A2 + A3)^{coskx - cosfcy)^ 



-(Al -I- ^2 - ^3)^ sin^ kx sin^ ky 



1/2 



,(47) 



which has been also obtained using displacement vari- 
ables in Ref. [T^ . Furthermore, from the Lagrangian 
equations, we find normalized symmetry-mode squared 



amplitude for the phonons of the square lattice with a 
monatomic basis, which represents how square plaque- 
ttes of four atoms distort when the lattice vibrates in 
each phonon mode. Corresponding expressions for gen- 
eral case, as well as for three special cases, are shown in 
Table U where Pi — l~ cos kx cos ky, (32 — ~ sin kx sin fcj,, 
= coskx — cos ky, Pi = {1 — cosfca;)(l — cosky), 
P5 = (l+cosfc^)(l+cosfcy), anda = {Ai-A2+A3)/{Ai + 
^2 — ^43). These expressions are also plotted within the 
first Brillouin zone of the square lattice in Figs. [HI [71 [51 
[HI and [TUl for a = 0, 0.1, 1, 10, and 00, respectively. The 
modes represented by the twelve panels in each figure are 
given in the caption for Fig. [6l 

We discuss overall features of these results. First, the 
squared amplitudes for the upper and lower branches 
show symmetries. The squared amplitude of the ei mode 
in the upper (lower) branch is identical to the squared 
amplitude of r mode in the lower (upper) branch (see 
general expressions in Table [I|. Similar relations exist 
between 62 and 63, between Sx and Sy, between tx and 
ty, between s+ and s_, and between i+ and Sec- 
ond, the squared amplitude is highly anisotropic. For 
a = 0.1 — 00, modes ei and 63 in the upper branch 
(equivalently, r and 62 in the lower branch) are dominant 
around k = (±7r,0), (0,±7r), and around the bond direc- 
tions from the zone center [panel (a)'s and (c)'s in Figs.[7| 
- 110], In contrast, mode 62 in the upper branch (equiva- 
lently, mode 63 in the lower branch) is dominant mainly 
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FIG. 6: Normalized symmetry-mode squared amplitude for 
phonons within the first Brillouin zone of the square lattice for 
a = 0, where a = {Ai - A2 + A3)/{Ai + A2 - A3). For the up- 
per branch, the twelve panels correspond to (a)|ei|^, (b)|e2p, 
(c)le3|^ (d)|rp, (e)ls.p, (f)l.,p, {g)\uf, (h)\ty\\ (i)l.+ |^ 
(j)|s_p, (k)|f-|-p, and (l)|t_p. For the lower branch, they cor- 
respond to (a)irp, (b)le3p, (c)le2p, (d)\er\\ (e)|s,p, (f)|s.p, 
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FIG. 7: Normalized symmetry-mode squared amplitude for 
the phonons for a — 0.1. Refer to Fig. |6]for the description 
of each panel. 
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FIG. 8: Normalized symmetry-mode squared amplitude for 
the phonons for a — 1. Refer to Fig. |6] for the description of 
each panel. 



in the region around zone-diagonal directions, which be- 
comes narrower as the parameter a increases [panel (b)'s 
in Figs. [7]- [To]. Third, the squared amplitude also de- 
pends strongly on the magnitude of wavevector k. Near 
the zone center, where fc w 0, tx, ty, t-^-, and i_ modes are 
dominant [panel (g)'s, (h)'s, (k)'s, and (l)'s in Figs. [S]- 
llOj . whereas near the zone corners, where k « (±7r, ±7r), 
Sx, Sy, s+, and s_ modes are dominant [panel (e)'s, (f)'s, 
(i)'s, and (j)'s in Figs. [HI -[TO]- Generally, the squared am- 
plitudes for Sx, Sy, s_(- and s_ at k are identical to those 



ig)\ty\'\ (h)|t.P, (i)|. 



P, (j)|.+ p, (k)|t_p, and (l)|t+|^ 



for tx, ty, t^, and i_ at (tt, tt) — k, respectively. Between 
/c w and k w (±7r, ±7r), ei, 62, 63, and r mode am- 
plitudes are prominent [panel (a) - (d)'s in Figs. [51-llOj. 
Fourth, the squared amplitudes depend on the moduli, 
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Mode 


General expression 


a = 


a 


= 1 


a = 00 




upper/lower branch 


upper/lower branch 


upper branch 


lower branch 


upper/lower branch 


|eiP 


1 ^ Sl+-Pl \ 

\ V ^2^" '-'Z / 




f3i 
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i(/3i±lAi) 


leap 


1 ^ /3?-aA? ^ 
« 1 Mi 1 9 ,02 i 
V V'32+" /^S / 


|(/3i±l/32l) 


/3i 
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i(/3iTl/33l) 


leap 


1 /3i-a;3i \ 
\ V P2+°- Pi J 




4/3i 


3S 


i(/3i±l/?3i) 




1 /3i+a,3i \ 

8 I ^1 + / 0-2. ^20-2 
\ V ' 2 ^ P3 / 


|(/3iTl/32l) 
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1/331) 




l^i fl-r "/Ss ^ 
8 \ ~r /fl2 I „2 fl2 1 
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^(/3i-/?3) 


^(/3i+/33) 


f [lTsign(/33)] 




Pi I t ^ 1/33 1 
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04 
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Wt - /33) 


f [l±sign(/33)] 
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i|[lTsign(/33)] 
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i^(/3i-/33) 


f [l±sign(/33)] 
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04 ( 1 ^ /32 1 
8 I ,/fl24_„2a2 I 
V V'^2^^ P3 / 


f [1 T sign(/32)] 






(34 
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fii I ^ P2 \ 

8 \ //32 1 2 i^2 1 


f [1 ± sign(/32)] 






(34 
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f [1 T sign(/32)l 






& 
8 




^ 1 ± "2 ^ 
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TABLE I: Normalized symmetry-mode squared amplitude for phonons within the first Brillouin zone of the square lattice. For 
± and =p, the upper sign corresponds to the upper branch, and the lower sign the lower branch. /3i = 1 — cos fcj, cos fcj, , P2 = 

— sinfca:Sinfcy, ^3 = cos fc^; — cos fcy , ^4, = (1 — cos — cos fcy), = (l+COsfca:)(l + COsfcj,), O = ( Ai — ^2 + ^3 ) / (^1 + ^2 " ^3 ) , 

and sign(a::)=x/|a;|. The amplitudes are normalized by |eip + |e2p + |e3p + |rp + Isvp + |syp + \tx\'^ + ~ 1. 



Ai, A2, and A3, through the parameter a, but are in- 
dependent of S. In particular, the general expressions 
bring additional symmetries for special cases of a = , 
1, and 00, which correspond to Ai + A3 — A2, A2 — A3, 
and Ai+ A2 ~ A3, respectively (Table|l|. Mode r in the 
upper and ei in the lower branches vanish completely for 
A2 — A3, which signifies the symmetry discussed below. 

We further comment on two aspects of the results, spe- 
cific to certain modes or values of a. First, as mentioned 
above, if a = 1, or A2 = A3, the upper phonon branch 
includes no rotational mode r, and the lower branch in- 
cludes no area-changing mode ei, which is interpreted 
in the following way. If the two shape-changing modes, 
62 and 63, have identical moduli, the lattice sustains 
isotropic phonon dispersion in the long wavelength limit, 
in which the lattice behaves like an isotropic continu- 
ous medium. Such medium would support longitudi- 
nal phonon modes in the upper branch and transverse 
phonon modes in the lower branch: the former rota- 
tionless and the latter locally area-preserving. For fi- 
nite wavelengths, the phonon dispersion is not isotropic, 
and the phonon modes are not longitudinal, nor trans- 
verse. However, the upper and lower branches of phonon 



modes remain locally rotationless and area-preserving re- 
spectively even for finite wavelength, if the two shape- 
changing modes have identical moduli. Panel (d)'s in 
Figs. [3 [SI and O show that amplitudes for these modes 
remain relatively small for a wide range of a. Sec- 
ond, the upper and lower branches of phonon modes 
at fc = (tt, 0) are longitudinal and transverse, as shown 
in Figs. [TIT a) and [TlTb). respectively. Therefore, the 
phonon at fc = (tt, 0) in the upper branch consists of ei 
and 63 modes, and that in the lower branch of 62 and 
r modes, which explains the large contribution of these 
modes near fc — (±7r, 0) and (0, ±7r) and around the bond 
direction from the zone center [panel (a)'s and (c)'s in 
Figs. [U- Hn]. For a ^ 0, or Ai + A3 ^ A2, the two 
phonon modes at fc = (tt, 0) shown in Fig. [TT] have the 
same net modulus and therefore are degenerate, which 
gives rise to the equal contribution of the four modes, ei, 
62, 63 and r at this fc point, as shown in panels (a) - (d) 
in Fig. E 

Understanding how different modes contribute to var- 
ious parts in the fc-space could be useful, for example, 
to gain insight in materials with electron-phonon cou- 
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FIG. 9: Normalized symmetry-mode squared amplitude for 
the phonons for a — 10. Refer to Fig. [6] for the description of 
each panel. 
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FIG. 10: Normalized symmetry-mode squared amplitude for 
the phonons for a — oo. Refer to Fig. [6] for the description of 
each panel. 
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FIG. 11: (a) Upper and (b) lower branch phonon modes 
at k = (7r,0) for the two-dimensional square lattice with a 
monatomic basis. 



pling, such as, manganites, phonon-mediated supercon- 
ductors, and materials near structural phase transition. 
The above analysis demonstrates the advantage of the 
approach based on symmetry mode for atomic scale de- 
scription of lattice dynamics. 



III. QUANTUM MECHANICAL FORMALISM 

A. One-dimensional lattice with a monatomic basis 

It is necessary to consider quantum mechanical aspects 
of lattice dynamics for phenomena such as low temper- 
ature specific heat, electron- phonon interaction, and po- 
larons. In this section, we extend the symmetry-based 
atomic scale description of lattice dynamics to the quan- 
tum mechanical formalism. Commutation relations be- 
tween the coordinate operators and their conjugate mo- 
mentum operators lie at the core of quantum mechanics, 
which we establish here for the symmetry modes. 

First, we consider the one-dimensional chain studied in 
Section II. A. The conjugate momenta for the two modes, 
Pe{i) and Pt{i) are 



dL 



M ., ^ 1 
— t{i) = — ■= 



{p,+i+Pi), (49) 



where pi represents the momentum of the atom at site i. 
From known commutation relations between momentum 
and displacement operators pi and Uj, we find the fol- 
lowing commutation relations between the operators for 
modes and their conjugate momenta with the same site 
index i, 

where a, 6 g {e, i}. The 1/2 factor is related to the num- 
ber of atoms in each motif. Unlike displacement vari- 
ables, the commutation relation between a mode at i and 
a conjugate momentum at j -I- 1 or i — 1 is not zero, since 
they share an atom, as shown below. 



[Pt(z),i(i±l)] = [Pe(*),i(* + 1)] = [Pt(z),e(i-1)] 



[PS). e{i ± 1)] = e{i + 1)] = [Pe{r),t{^ - 1)] 



h 

4i' 



4i 



The commutation relations between the momentum and 
the mode, defined at sites further than the nearest neigh- 
bors, vanish. 

The above relations are also established graphically. 
For example, [Pe(«),i(* + 1)] is found from the drawing 
in Fig. [T21 where Pe{i) and t{i + 1) are represented with 
arrows. We treat the arrows as unit vectors, and find the 
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i i+l i+2 
- • •-» 

Pe(i) t(i+l) 

FIG. 12; Graph to find commutation relation, [Pe(i), t(i + 1)]- 



sum of scalar products of unit vectors defined at the same 
sites, which after being multiphed by fi/(2^i), lead to the 
commutation relation. From the graphical rule and the 
symmetry of the modes, the following commutations are 
obtained, where a,b € {e,t}. 

[Pa{t),b{j)] = [Pb{j),a{t)], (50) 
[Pe{i),m = -[Pt{^),e{j)]. (51) 

The commutation relations in reciprocal space are calcu- 
lated from the relations. 



[F,(fc),6(fc')] = Sk',-kJ2iP-i' = O)'^(J)]' 



,ikj 



(52) 



which are shown in Table |TT1 The results found here 
are applicable, for example, for the study of quantum 
mechanical dynamics of non-linear excitations mentioned 
in Section II. A. 



modes are as follows. 

Peiii) = ^ei(i) ^ 



Pr& 



PsA^)-- 

PtJ) 
Pty& 



M 



e2& 



M 



r{i) 



4 
M 



4 
M. 

M. 



1- f-p5-pH 



8\/2 
1 / 



y 



8^/2 



-Pr 



pi 



i+Ol 



8\/2 



-P7 



1 / 



y 
"P- 



-PUii 



V- 



-PX+ii 



y 



- pH 



■ pH 



y 
p" 



y 



8V2 



V 



y 
p- 



-Pu 



Sy{i) = 



8l^? 
8 
8 
8 



' Pi+io 



y 
p- 



- P- 



-^4+ 10 



y 
p" 



- P- 



■ pi 



' pi 
^i + 11 



-Pr+11 



pi 

^i+ll 



-PUii 



- pi 



From the fundamental commutation relations for dis- 
placement operators and momentum operators. 





iPe(fc) 




e(-fc) 


1(1 — COS k) 


— 1 sin k 


t(-fc) 


1 sin k 


|(1 + cos k) 



[p^.,U^^[pi,ul] = -^77 



TABLE II: Commutation relation, [^Pa(fc), 6(— fc)], between 
symmetry modes and their conjugate momenta for the one- 
dimensional chain in reciprocal space. 



B. Two-dimensional square lattice with a 
monatomic basis 

Quantum mechanical nature of lattice is also impor- 
tant for two or three-dimensional lattices, for example, 
near the structural phase transitions. In this subsection, 
we find quantum mechanical commutation relations for 
the symmetry modes and their conjugate momenta for 
the square lattice studied in Section II. B. 

Conjugate momenta for the atomic scale symmetry 



[p5, ui 



hi' 



the commutation relations between modes and their con- 
jugate momenta are calculated in a straight forward way. 

However, it is more convenient to use the graphical 
method, explained for the one-dimensional chain in the 
previous subsection. The above fundamental commuta- 
tion relations for i — j have the form of 

x-x^y-y = 1, 
x-y^y-x = 0, 

except for the factor h/i, where i, and y represent unit 
vectors, not operators. Therefore, the commutation rela- 
tion [Pa(«), ^(j)]) where a and h represent the eight atomic 
scale modes, is found from the drawings of a and h modes 
on the square lattice. The sum of the scalar products of 
the unit vectors at the sites shared by the two modes, 
multiplied by h/{A^i), gives the commutation of the two 
operators. (The multiplication factor after h/i is associ- 
ated with the number of atoms in the motif for the lattice 
with a monatomic basis, that is 4 for the square lattice 
and 2 for the chain.) For example, [Pei(i),e2(i -I- 11)] is 
found from Fig. [T31 as follows 



[P,i(*),e2(z-l-ll)]-^(-l). 



(53) 
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-^Pe3{k) 




i r) /IN 


in / j„ \ 




nPtyW 


ei(-fc) 






~^kx'^^ky 





i{l-Cfc^)Sfc„ 


i(l-Cfc„)Sfc^ 




i{l+C^^)Sky 


4 


4 


4 




4\/2 


-4\/2 


-4V2 


e2(-fc) 




^ h X ^ ky 





^ re 37 1 ^ Ky 


i(l — Ck„)Sk^ 


Z f X t. T ) 1 < 

V ^ J ^ tuy 


i ( 1-l-C' I S u^. 

\ 1 ^ KH^ / ^ Ky 


i(l-\-Cii.^,)Sh.^ 

\ ' ^ Ky J ^ KX 


4 


4 


4 


A-J2 


4\/2 


-4s/2 


-4V2 






U 


^~^kx ^ky 




^{^-Ckx)Sky 




i(i+Cky)Skx 




4 


4 


4 




-4v'2 


-4\/2 


4^/2 


T\—k) 





~^'kx '^^'ky 


^ kx ^ ky 


'^ — ^kx ^ky 


i(l-Cky)Skx 


)Sky 


i(l+Cfcx)Sfcy 


^{'^ + Cky)Skx 


4 


4 


4 




4V2 


4v^ 


-4V2 


Si:( — fe) 


i(l-Cfcx)Sfcj 








(1-Cfc^){l-Cfc^) 





— '^kx ^ky 





-4\/2 


-4^ 


-4^ 


4^ 


4 


4 




«(l-Cfcj,)Si.j. 




i(l-Cfc„)Sfc, 


i{l-Ckx)Sky 


n 
u 


{1-Cfc,)(l-Cfc„) 


n 
U 


-^kx^ky 


-4v^ 


-4v^ 




-4V2 


4 


4 










i{l + Ckx)Sky 


'^kx '-'ky 





(1+Cfcx)(l+Cfcy) 





4v^ 


4v^ 


4v^ 


-A-J2 


4 


4 




s(l + Cfcj;)St.„ 






i(\ + Cky)Skx 





— Skx ^ky 





(l + CkxK-^ + Cky) 


4^2 


4^2 


-4^2 


4,V2 


4 


4 



TABLE III: Commutation relation, [jPa{k),b{—k)], between symmetry modes and their conjugate momenta for the two- 
dimensional square lattice in reciprocal space. Ckx, Cky, Skx, and Sky represent cosfcj,, cos ky, sinfca;, and sinfc^, respectively. 



o o 

^ \ 

FIG. 13: Graph to find commutation relation, [Pei(j),e2(i + 
11)]. 



Presented graphical method is also useful to find the 
following symmetry-related properties of the commuta- 
tion relations, where a and b represent any of the eight 
modes, and even and odd represent the modes with even 
symmetry under point reflection, namely, ei, 62, 63, r, 
and the modes with odd symmetry, namely, Sx, Sy, t^, 
ty, respectively. 



[Pa(^)Mi)] 

\Peven ('^) ; 

even'{j)] 

[Peven{i),odd{j)] 
[Podd&,0dd'{j)] 

The commutation relations in reciprocal space are found 



[Peven'ii),even{j)], 
-[Podd{i),even{j)], 
[Podd' {i),odd{j)]. 



from the relation 

[Pa{k),b(k')] = 5- _-^[P,(r= 0), 6(J)]e'^'-? (54) 



which are provided in Table IIIII 

IV. SUMMARY 

In this article, we have presented mode-based atomic 
scale description of lattice dynamics. It is found that not 
only the potential energy but also the kinetic energy is de- 
scribed in terms of the atomic scale modes, for which the 
inclusion of the rigid modes is essential. This approach 
has been demonstrated for the one-dimensional chain and 
the two-dimensional square lattice with a monatomic ba- 
sis. The comparison with a continuum model has shown 
that our approach is suitable for multiscale description 
of lattice dynamics. Phonon analysis in terms of sym- 
metry modes have revealed the momentum dependence 
associated with underlying symmetry of the modes. The 
approach has been extended to quantum mechanics, and 
the commutation relations have been obtained. We ex- 
pect that this method would be useful in describing 
atomic scale lattice dynamics in systems with strong an- 
harmonicity and complex energy landscape, which can 
be compared with the results from time-resolved x-ray 
experiments. 
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